Transverse laser modes in Bose-Einstein condensates 
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We examine the Bose-Einstein-condensate (BEC) equivalent of transverse aspects of laser res- 
onators. We model numerically repeated focussing of a 2-dimensional BEC, which could be achieved 
in practice by a series of far off-resonant light pulses. We show for a range of non-linear coefficients 
that such a series of light pulses traps the BEC. We also model a combination of repeated focussing 
and loss, which could be achieved with a series of light pulses, some far off-resonant, some res- 
onant. In analogy to light repeatedly being focussed and locally absorbed by passing through a 
laser resonator, a small proportion (in our model up to 10%) of the BEC is shaped into Hermite- 
Gaussian-like modes; the remainder of the BEC is lost. This happens irrespective of the presence of 
a harmonic trap. We show that astigmatic focussing of the resulting Hermite-Gaussian-like modes, 
which can be achieved by a pair of off-resonant light pulses, results in a specific number of vortices. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) are the matter- 
wave equivalent of lasers p): all the atoms in a BEC 
are in the same quantum state, just like all the photons 
in a laser beam. While in the creation of BECs this prop- 
erty is usually ensured by the fundamental properties of 
very cold bosons, in a laser it is due to the longitudinal 
properties of the laser resonator Just like the num- 
ber of photons can be coherently amplified in a laser, 
the number of atoms can be coherently amplified in a 
BEC HH0. 

Laser resonators also have an effect on the transverse 
structure of laser light. The origin of this transverse 
structure has a huge literature devoted to it (for an 
overview see, for example, references 0,0); and a good 
understanding and control has proved interesting (see, 
for example, the discovery of fractal modes Q) as well 
as useful. For example, the standard families of laser 
modes, Hermite- Gaussian (HG) and Laguerre-Gaussian 
(LG) modes (figure 0, have elegant mathematical prop- 
erties (for example, they are structurally stable, i.e. they 
do not change shape on propagation, and form complete 
orthogonal sets of modes), and LG modes, because of 
their central vortex of arbitrary topological charge, are 
one of the most popular choices for optical experiments 
with optical vortices. Most importantly, however, certain 
'stable' resonators keep the light in the proximity of the 
optical axis - transversally, they trap the light. 

We examine here - in the limit of the Gross-Pitaevskii 
equation - a process in 2-dimensional BECs that is anal- 
ogous to the formation of transverse laser modes: peri- 
odic focussing and spatial filtering. We believe this could 
be achieved experimentally with a periodic series of light 
pulses. We examine this scenario for BECs in a quadratic 
trap and untrapped BECs. We demonstrate the emer- 
gence of HG-like modes, and convert these into LG-likc 
modes using the equivalent of the Gouy phase shift. Our 
main concern, however, is the question if a series of light 
pulses could act as a novel type of trap for BECs. 
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FIG. 1: Hermite- Gaussian (HG) and Laguerre-Gaussian (LG) 
modes. Shown here are the intensity cross-sections of two 
HG modes, HG2,o and HG2,i, and the intensity and phase 
(grayscale representation) cross-sections of one LG mode, 
LG4,o, which possesses a central charge-4 vortex. 



II. EIGENMODE SELECTION IN LASER 
RESONATORS 

Figure Eh sketches one round trip of a light pulse (or 
alternatively a 'slice' of a cw light beam) through a stan- 
dard laser resonator. Most importantly, the light pulse 
is reflected back by the two mirrors facing each other, 
so that after one round trip, which involves reflecting 
off each mirror once, the light pulse is back in the same 
plane where it started. One round trip is much more 
complicated than just some changes in propagation di- 
rection: as the mirrors are curved, each mirror reflection 
also focusses the light. In addition, localised loss occurs 
in one or several planes, for example through apertures 
or specks of dust on mirrors. Some light beams - the res- 
onator's eigenmodes - are unchanged (apart from a uni- 
form change in intensity or phase) after one round trip 
through a resonator. In the case of geometrically stable 
resonators, which we consider here, the eigenmodes can 
be divided into families (like the Hermite-Gaussian (HG) 
modes and the Laguerre-Gaussian (LG) modes), whereby 
each family forms a complete orthogonal set, i.e. every 
light beam can, for example, be described as a super- 
position of HG modes. As already mentioned, the HG 
and LG modes can take on very recognisable shapes and 
have interesting features like high-charge vortices (see fig- 
ure nj. 
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FIG. 2: Laser resonator (a) and its unfolded lens-series equiv- 
alent (b). A light beam passing through the lens-series equiv- 
alent is periodically being focussed, just like light in the res- 
onator is being focussed every time it is reflected off one of 
the resonator's mirrors. A series of 'lens pulses' - short, off- 
resonant light pulses with a specific intensity profile - can 
repeatedly focus a BEC in a similar way (c). Note the role 
reversal: in a and b, a light pulse is being focussed (by mirrors 
and lenses, respectively); in c light pulses focus the BEC. 



Experimentally, a pure resonator eigenmode is cre- 
ated in a laser by making that eigenmode the lowest- 
loss eigenmode, i.e. by ensuring the round-trip loss of 
all other eigenmodes is higher than that of the desired 
eigenmode. This can be done, for example, by insert- 
ing absorbing cross-wires into the resonator at positions 
where the desired eigenmode is darker than the compet- 
ing eigenmodes. For example, an absorbing wire in the 
central horizontal dark line of the HG2,i mode shown in 
figure^would not significantly increase that eigenmode's 
round-trip loss, but it would increase that of the HG2,o 



J 



mode. The fraction of the power in the undesired eigen- 
modes then falls of exponentially with round trip number, 
after a few round trips leaving only the pure lowest-loss 
eigenmode. 



III. CREATION OF 
HERMITE-GAUSSIAN-LIKE MODES IN BECS 



As far as transverse eigenmode selection is concerned, a 
resonator can be seen as a series of lenses separated by the 
length of the resonator (figure [2b) ■ If the focal lengths 
of the lenses match those of the resonator's mirrors, and 
if the important localised loss is the same, then an in- 
finitely long lens-series equivalent of a resonator has the 
same eigenmodes 0. Such an infinitely long lens-series 
equivalent of a resonator represents infinitely many round 
trips. On the other hand, eigenmode selection mostly 
happens during the first few round trips through the res- 
onator, which can be represented by a few lenses. We 
examine here the BEC analog of such a series of lenses. 

The BEC analog of an optical resonator's lens-series 
equivalent can consist of repetitions of lens pulses (fo- 
cal time /) and absorption events (figure 0:; absorption 
events are not shown). The delay between two repeti- 
tions, L, corresponds to the length of the resonator. We 
refer to such a configuration as a light-pulse resonator. 
We restrict ourselves to light-pulse resonators in which 
the absorption events always follows immediately after 
the lens pulse. 

All the pictures of BECs in light-pulse resonators 
shown in this paper (with the exception of two frames 
in figure Q) are snapshots taken half-time between two 
lens pulses, i.e. a time L/2 after the previous lens pulse. 
We refer to the sequence of events acting on the BEC 
between L/2 after the previous lens pulse and L/2 af- 
ter the next lens pulse - time evolution over a time L/2, 
interaction with a lens pulse, absorption event, and an- 
other time evolution through L/2 - as one period of the 
resonator. 

We simulate time evolution according to the time- 
dependent Gross-Pitaevskii equation [9| , which we write 
in the form 



1 dt 



fir i 

- — V 2 + -mu?r 2 + V e (r,t) + g\i/>(r,t)\' 
2m 2 



li>(T,t), 



(1) 



r 



where m is the mass of each atom, u> is the trap frequency, 
V e is an external potential, and g is the non-linear coef- 
ficient, which is given by 



g = AttN 



d m 



(2) 



d is the thickness of the BEC in the third dimension (the 
z direction), which needs to be very small for the BEC 
to behave 2-dimensionally. This is usually ensured by 
tight z confinement by a harmonic trap that satisfies the 
condition fi < huj z , where /i is the chemical potential. 
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FIG. 3: Focussing of a BEC by a lens pulse - a far off-resonant 
light pulse with a quadratic intensity profile that phase-shifts 
the BEC locally by r 2 /(2f). The graph shows the numeri- 
cally calculated relationship between the 'focussing time' pa- 
rameter, /, and the actual focussing time, tf, after which a 
BEC that was initially in the trap ground state is found to be 
phase- flat. Analogy with optics suggests that the BEC comes 
to a focus after a time /, i.e. tf = f; however, to our surprise 
we find that this is only the case for small values of /. Our 
simulations suggest that the relationship shown in the graph 
holds for values of g in the range < g < 100. 

Note that g depends not only on m and a, the s-wave 
scattering length, but also on the number of atoms, N, 
in the BEC. As the number of particles is contained in the 
parameter g, which is very convenient for the purposes 
of this paper, the density is normalised according to 

J \^(y,t)\ 2 dv = l. (3) 

In our computer model we represent the BEC wave func- 
tion on a square area of 16 x 16 dimensionless units by 
a 128 x 128 array of double-precision complex numbers. 
The dimensionless variables for position and time, r' and 
t' , are related to the corresponding variables with dimen- 
sions, r and t, through the equations 

*={—) r (4) 

and 

t' = ut. (5) 

We solve the Gross-Pitaevskii equation with a fourth- 
order Runge-Kutta method (a detailed description can 
be found in ref. [To|). 

We model the effect of a lens pulse by locally alter- 
ing the phase of the wave function in the same way in 



which a lens (or a hologram of a lens) would alter the 
phase of passing light [Tj. In principle (although it has 
not been demonstrated for the particular phase profile 
of a lens), this can be achieved experimentally using an 
off-resonant laser pulse with an intensity profile propor- 
tional to the desired local phase change, whose electric 
field then - through the AC Stark effect - alters the lo- 
cal phase of the BEC. This 'phase imprinting' method 
[13. Il3| is analogous to optical phase holography. Inter- 
estingly, the time after which the BEC comes to a focus 
after interaction with such a lens pulse is not what one 
might expect it to be from the optical analogy (see figure 
[3] and appendix). 

At this point, it is perhaps worth to briefly discuss the 
BEC-light analogy in slightly more detail. The time evo- 
lution of an interaction-free, un-trapped, 2-dimcnsional, 
BEC according to the Gross-Pitaevskii equation is for- 
mally equivalent to spatial propagation of a light beam 
in the paraxial approximation. This is, of course, the 
basis of our specific efforts to reproduce light-like be- 
haviour in BECs. However, the analogous light beam 
is of sub-wavelength size, and the paraxial approxima- 
tion describes the actual propagation of such beams very 
badly (it does, for example, ignore evanescent waves). It 
is therefore not surprising that BECs - even in the most 
light-like case of no interactions and no trap - often be- 
have quite differently from light. 

We model local absorption as a straightforward reduc- 
tion in amplitude; in particular we neglect any interac- 
tion between the 'absorbed' and left-over parts of the 
BEC. This is an idealisation of a situation which, wc 
believe, can be approximated experimentally in differ- 
ent ways, for example by utilising a two-photon Raman 
process to couple out (move into an untrapped state) 
the part of the BEC which lies in the overlap region of 
two laser pulses [l4| . To achieve the necessary flexibility, 
our absorption pattern consists of up to two absorbing 
cross-hairs (the intensity cross-section is that of a narrow 
Gaussian, in our simulation often a single pixel wide; for 
a typical BEC size of 50/im, this corresponds to lines of 
light about 500nm wide) , as well as an absorbing rectan- 
gular boundary with a Gaussian edge profile. 

Figurc0|shows the BEC after the first few light pulses, 
starting with the trap ground state, through a light-pulse 
resonator of length L = 0.2 and with a focal time / = 0.5 
(both in dimensionless time units). We use this as our 
'standard' light-pulse resonator. It will be seen that this 
behaves in many ways analogous to a stable optical res- 
onator. Due to the losses, the number of particles changes 
during every period of the light-pulse resonator, which in 
turn changes the size of the self-interaction potential in 
the Gross-Pitaevskii equation, g\ip\ 2 (provided g 7^ 0). 
As we keep the integral over the density normalised in 
our model, the reduction in the number of particles and 
any corresponding change in the interaction term size is 
represented as a decrease in g. 

As long as the value of g - and with it the potential - 
changes, no stationary eigenmode can form. This means 
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FIG. 4: BEC after interaction with p lens pulses and absorp- 
tion events of a light-pulse resonator (the value of p for each 
frame is shown in the top-left corner), starting with the har- 
monic trap ground state (p = 0) , calculated for g = 100. The 
pulses of focal time / = 0.5 are separated by L = 0.2 (in 
dimensionless time units); the density patterns shown here 
are calculated a time L/2 after the actual pulse, i.e. half-time 
between the previous and the next pulse. The curve at the 
bottom shows the fraction of the initial number of particles, 
No , that remains in the BEC after p periods of the light-pulse 
resonator. To demonstrate the loss, the density is not nor- 
malised in the frames in this figure; the full scale of grays 
represents in the frames for p = 1 to 4 densities between 
and 0.2, in the remaining frames between and 0.05. 



that, for g ^ 0, no stationary eigenmode can form be- 
fore the loss has become zero. In practice, the number 
of particles keeps diminishing exponentially, so zero loss 
can never be reached with any atoms still left. However, 
while atoms are still left, the loss can become very small. 
Then the BEC becomes very similar to the eigenmode 
corresponding to the eigenmode for the same value of g 
in the hypothetical limit of no loss. 

In order to calculate eigenmodes for any given value of 
g in the no-loss limit, we re-normalise the integral over 
the density after each loss pulse without making a cor- 
responding adjustment to g. This would correspond to 
either coherently replenishing the BEC after each loss 
pulse such that the number of atoms remains constant, 
or to tuning the value of g to 'simulate' this - neither 
of these are easily realisable experimentally. Neverthe- 
less, the procedure gives stationary eigenmodes that cor- 
respond to the no-loss limit. 

After a number of periods of the light-pulse resonator 
the eigenmode has formed, that is the BEC is not signifi- 
cantly changing any more over subsequent periods (apart 
from a uniform factor to the wave function). As a cri- 
terion for reaching this state we demand that the loss 



FIG. 5: Density patterns of some eigenmodes of the stan- 
dard light-pulse resonator with a non-linearity of g = 30. By 
varying the cross-hair position, different patterns - charac- 
terised here by their number of nodal lines in the horizontal 
and vertical direction, n and m, respectively - emerge, which 
strongly (but not exactly) resemble optical Hermite-Gaussian 
modes. These pictures show the BEC half-way between two 
lens pulses. A possible position of the cross-hair (white dashed 
lines) is shown in the two rightmost frames. The pictures show 
only the central 6.25 x 6.25 units of the modelled area. 



due to the aperture varied by less than 0.5%; this hap- 
pened typically after about 40 iterations (light-pulse se- 
quences), which takes several minutes to simulate on our 
desktop computer. Note that there appears to be an up- 
per limit to the non-linear coefficient, g m a.x, above which 
the BEC does not converge to a stable eigenmode. The 
exact value of <7 ma x depends on the parameters of the 
light-pulse resonator. We determined this value for a 
number of different resonators and obtained values be- 
tween 80 < <7max ^ 110. 

Once the eigenmode has formed, the absorption due to 
each absorption event can be very low; we have achieved 
losses below 0.2% per absorption event. However, before 
the eigenmode has formed the losses can be very high; 
just how high depends on the starting conditions and 
the particular choice of shape of the loss pulse. Note that 
there is a trade-off between round-trip loss and number 
of round trips it takes to form eigenmode: thinner cross 
wires mean lower loss per round trip, but require more 
round trips. For well-chosen parameters, of the order of 
10% of the BEC are left by the time the BEC begins to 
resemble the stationary eigenmode. 

Figure[S]shows some of the stable patterns that emerge 
for a modest non-linearity (g = 30, corresponding to ap- 
proximately 10 3 atoms of s7 Rb in a 20tt Hz trap) by vary- 
ing the localised absorption pattern - a 100% absorptive 
cross-hair of 1 pixel width placed in crossing node lines 
(see figure |5J. Like all other pictures in this paper (un- 
less otherwise stated), this figure shows only the central 
6.25 x 6.25 units of the modelled area for clarity, and 
the grey values represent linear densities, ranging from 
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FIG. 6: Density (top) and corresponding phase patterns (be- 
low) for various values of the non- linear coefficient, g. In all 
cases, the shape of the absorption region is a narrow vertical 
straight line. 



(black) to 0.2 (white). The patterns look very similar to 
Hermite-Gaussian modes, more so for small values of the 
non-linear parameter g, less so for large values (figure|BJ). 
In all cases, the eigenmodes have the tell-tale ir phase 
difference between adjacent peaks, just like the optical 
HG modes. 

Repeating these simulations for a switched-off har- 
monic trap shows that the results are virtually indistin- 
guishable from those obtained with the trap switched on. 
This is not surprising as the time scale of effects due to 
the trap (1 dimensionless unit) is several times larger 
than the time scale of effects due to the lens pulses (the 
focal time). It suggests that a series of lens pulses - a 
light-pulse resonator without absorption events - could 
act as a novel type of trap. Indeed, in the absence of ab- 
sorption events, and to within the accuracy of the Gross- 
Pitaevskii equation, a series of light pulses acts as a loss- 
less trap. It is worth noting that, much like a stable laser 
resonator 'traps' not only light waves but also light rays, 
a light-pulse resonator could perhaps not only trap atoms 
in a BEC but also hotter or Fermionic atoms. 

If the series of light pulses, as well as any other po- 
tential (in particular the harmonic trap), is switched off 
after a stationary eigenmode has formed, the eigenmode 
evolves freely. Figure investigates the size change of 
the eigenmode during this free evolution; the shape re- 
mains almost unchanged. This property is in analogy to 
structural stability of the eigenmodes of stable optical 
resonators, which can be understood in terms of imaging 
of every plane into every other plane [T^ . 

All the results shown so far were calculated for our 
standard resonator with a time between pulses - cor- 
responding to the resonator length - of L = 0.2, and 
a lens-pulse focal time of / = 0.5. Figure [H] shows 
n = 1, m = eigenmodes from different light-pulse res- 
onators. In the optical analogy, altering resonator pa- 
rameters while keeping the resonator geometrically sta- 



FIG. 7: Free time evolution of an eigenmode, calculated for 
g — 0. If - after the eigenmode has formed - the series of 
light pulses is switched off at t = (the last few pulses arrive 
at t = — 2L, —L, 0), the BEC approximately retains its shape 
but changes its size. The graph shows the distance from the 
centre of the maximum of the density, which indicates the rel- 
ative size; the points do not lie on a smooth curve due to nu- 
merical errors. In analogy with optical higher-order Hermite- 
Gaussian modes, the BEC passes through a 'focus', where it 
is smallest. On either side of the focus, its size changes ap- 
proximately symmetrically. The data points for t < were 
calculated using 'backward propagation' in time. The graph 
is calculated for L — 0.2 (in adimensional time units); the ver- 
tical range then represents 3 adimensional length units. To 
avoid numerical artefacts in obtaining this graph, the BEC 
was simulated over an area larger than the standard 16 x 16 
non-dimensional units. 
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FIG. 8: Density patterns for the output of different configu- 
rations of resonator, calculated for g = 0. / represents the 
focal time of the lens pulse, L the time between subsequent 
pulses (both in dimensionless units). For the particular choice 
of loss-pulse shape, the loss per period of the light-pulse res- 
onator for these eigenmodes is, from top left to bottom right, 
1.8%, 1.3%, 5.9%, 4.4%, 10.2%, and 9.7%. 
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FIG. 9: Modelled density (top) and phase (bottom; grey- 
level representation) after a pair of cylindrical-lens pulses 
has passed through eigenmodes with, from left to right, 
n = 1,2,4, respectively, and m = 0, in a BEC with g — 50. 
The patterns were obtained with cylindrical-lens pulse pairs 
with, from left to right, focal time / = 0.23 and separation 
s = 0.35, / = 0.23 and s = 0.3, and / = 0.1 and s = 0.11, 
respectively. The phase patterns show the presence of re- 
spectively I = 1,2, 4 vortices, situated at positions where all 
possible phases meet (marked by arrows). 

ble (that is, satisfying L < 4/) results in eigenmodes of 
a different size; the analogous effect in BEC light-pulse 
resonators can be seen in figure |S] 



IV. CONVERSION INTO 
LAGUERRE-GAUSSIAN-LIKE MODES 

To examine the similarities between the eigenmodes of 
light-pulse resonators and optical HG modes beyond the 
level of appearance, we investigate here the analogy of 
a conversion of optical HG modes, which have a rect- 
angular symmetry, into circularly symmetric Laguerre- 
Gaussian (LG) modes under cylindrical focussing [lfij . 
This conversion is due to a fairly subtle effect: the Gouy 
phase shift, which occurs when HG modes pass through 
a focus. It is interesting also because LG modes can con- 
tain a higher-charge vortex at their centre. 

In the optical case, the cylindrical focussing that con- 
verts HG modes into LG modes is usually performed with 
a pair of cylindrical lenses. In analogy, we use here a pair 
of cylindrical-lens pulses, i.e. far off-resonant laser pulses 
with an intensity profile that varies quadratically in one 
direction and which is constant in the perpendicular di- 
rection. FigureElshows some results that were calculated 
for g = 50. The density is certainly not shaped like a per- 
fect ring; instead, most of the density is contained in a 
distorted ring. Also instead of containing a charge-^ vor- 
tex at the centre, the wave function contains I vortices 
of charge 1. The conversion is therefore not perfect, but 
there are some indications that it is happening. In other 
simulations (not shown) we confirmed that this is also 



the case for different choices of g (including g — 0), and 
that, in the appropriate cases, the shape of the BEC also 
shows signs of the multiple concentric rings characteristic 
of LG modes with a radial mode index p > 0. 

In the optical case, the focal lengths of the cylindrical 
lenses and their positions need to be carefully matched to 
the parameters of the incoming HG mode to achieve good 
conversion into a LG mode. We tried to similarly opti- 
mise the parameters of the cylindrical-lens pulses. How- 
ever, the similarity to a pure LG mode of the resulting 
density was always limited, which is not surprising as the 
BEC was initially not in a pure HG mode and as the evo- 
lution of the BEC is different from that of light. Finally, 
optimisation of the pulse parameters was not made easier 
by the complications in the focussing characteristics of a 
BEC mentioned earlier (figure |3J. 



V. CONCLUSIONS 

The emergence of Hermite- Gaussian-like modes in a 
BEC is a very graphical further illustration of the sim- 
ilarity between BECs and laser light. However, any ex- 
perimental verification of our scheme would have to take 
into account the high losses associated with the absorp- 
tion events. 

Alternatively, the absorption events could be replaced 
by far off-resonant pulses that locally disturb the BEC 
in places where the local density of the desired mode is 
zero, thereby effectively transferring a small part of the 
BEC in an undesired mode into a variety of other modes, 
amongst them the desired mode. No simulations of this 
scheme have been performed to date, so the magnitude 
of any loss saving is unknown. 

Our results indicate that a series of lens pulses acts as 
a novel trap. Although analogy with optical resonators, 
in agreement with our approximate description with the 
Gross-Pitaevskii equation, would suggest that such a trap 
is loss-less, a more detailed microscopic description of 
the BEC predicts finite losses during the interaction with 
the lens pulses and when the density gets hi gh d ue to 
mechanisms such as 3-body recombination [171 

As described in this paper, a light-pulse resonator con- 
fines a BEC in two dimensions. It should be possible to 
achieve confinement in all three dimensions by a combi- 
nation of two orthogonal light-pulse resonators or similar 
arrangements. Alternatively, the rapid fall-off of the in- 
tensity away from the focal plane due to extremely tight 
focussing of the lens pulses, which leads to an approxi- 
mately quadratic Stark potential along the propagation 
axis, could perhaps be used to extend the confinement to 
all three dimensions. 

We are currently investigating the BEC analogy to 
holographic light shaping. In this way, a BEC can be 
shaped in two and three dimensions. This could, for ex- 
ample, be useful in quantum computation with arrays of 
trapped atoms 19] . The analogy can even be extended 
to the extremely versatile volume holograms, whereby 
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the BEC analog of a controllable volume phase hologram 
- a far off-resonant light pulse with a controllable time- 
varying intensity - might even be easier to realise than 
the original. 
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APPENDIX A: 



ANALYTICAL FOCUSSING 
RESULTS 



Let us briefly discuss the focussing dynamics in two 
dimensions. For low densities or small scattering lengths 
we can assume a non-interacting gas, i.e. g = 0. Using 
the propagator for the now linear Schrodinger equation 
with a harmonic external potential we obtain the full 
dynamics of the condensate, 

. / r 2 . ,\ l + icot(wi) 

\Er(r, t) = exp -i — cot(u)t) ; ; — 7 — —r-^ 

y ' ; V 4 /1+ i(cot(ujt) - 1//) 



x exp 



4shU(wt) (1 + i(cot{ut) - 1//)) 



where we have used the initial state 



*(r, t = 0) = exp 



1 2 ■ 1 2 
— r — 1 — r 
4 4/ 



(Al) 



(A2) 



This initial state is a ground state whose phase has been 
altered by a lens pulse of focussing time /. All lengths 
are in units of (h/2muj) 1 / 2 and the energy is in units of 
Tilj. From the minimum of the width of the Gaussian 
\^(r, t)\ 2 we obtain the actual focussing time 



1 1 

ui2 



arctan(2/). 



(A3) 



If / 3> 1 we get tf = 7r/4. In the opposite limit, / 1, 
tf w /, in accordance with figure |31 

We can also calculate an analytical expression for the 
focussing time when the interactions between the atoms 
are not negligible and the inequality fx/fiui ^> 1 (where 
fi is the chemical potential) is satisfied - the Thomas- 
Fermi regime. During the focussing the condensate is 
structurally stable, i.e. the shape of the cloud is preserved 



and only the time and length scales change during the 
evolution. We rewrite the Gross-Pitaevskii equation in 
the hydrodynamic form 



^+V-(pv)=0 
dv _/l 



(A4) 



dt 



= -V 



2 2 

-rauj r 



h 2 v ; 



gp 



2m ^fp 



(A5) 



where p(r, t) is the density with the velocity v defined as 
v = (ft/m)V9(r,t) and V6(r,i) is the gradient of the 
condensate phase. In the Thomas-Fermi approximation 
we neglect the last term on the right-hand side of Eq. 
±51. The density will now take the form 



p(r,t) 



^ [ »-l mu2{ W) n (A6) 



where p, is the chemical potential at t = 0. If we choose 
the phase to be quadratic in r, 



0(r,t) = —r 2 q(t) 



(A7) 



and insert this into Eqs. (|A4|) and (|A5|) together with 
the parabolic density p(r, t), we see that the dynamics is 
completely described by the parameter b(t) if and only if 
b(t) fulfills 



1 



(A8) 



with q(t) = b/b and the initial conditions 6(0) = 1 and 
6(0) = —1//. Equation (|A8|) can be solved exactly to 
give 



bit) = A/l + ^-|V^ + ^ Sin ( 2£ + arCtan( 27 } 



(A9) 



The focussing time, defined by 6 = 0, is consequently 



tf = — - arctan(2/) 



(A10) 



As can be seen from Eqs. (|A3() and (|A10|) the focussing 
time is the same both for the noninteracting case and 
in the Thomas- Fermi limit. This is a coincidence and 
holds only in the self-similar situation and in 2D since 
the radial excitation frequencies are given by = 2nu 
for the ideal gas and O = u>\/2n 2 + 2nm + 2n + m [2£j 
for the Tomas-Fermi cloud. For n = 1 corresponding to 
the lowest breathing mode, i.e. focusing, the frequencies 
are the same. 
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